> 
in 

00 

d 

(N 



X: 



The Equivalence of the Lagrangian-Averaged Navier-Stokes-a Model and the 
Rational LES model in Two Dimensions 

Balasubramanya T. Nadiga-^O and Freddy Bouchet^ l^ 

^^LANL, Los Alamos, NM-87545 
^>ENS-Lyon, CNRS, Lyon, France 

O 

^ ■ (Dated: 9 October 2012) 

O ■ In the Large Eddy Simulation (LES) framework for modeling a turbulent 
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flow, when the large scale velocity field is defined by low-pass filtering the 
^ ' full velocity field, a Taylor series expansion of the full velocity field in terms 

• of the large scale velocity field leads (at the leading order) to the nonlinear 



gradient model for the subfilter stresses. Motivated by the fact that while 



^3 

c/3 ■ the nonlinear gradient model shows excellent a priori agreement in resolved 

simulations, the use of this model by itself is problematic, we consider two 
jS^ . models that are related, but better behaved: The Rational LES model that 



uses a sub-diagonal Fade approximation instead of a Taylor series expansion 
and the Lagrangian Averaged Navier-Stokes-a model that uses a regulariza- 
tion approach to modeling turbulence. In this article, we show that these two 
latter models are identical in two dimensions. 
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I. INTRODUCTION 



In a turbulent flow, it is usually the case that energy is predominantly contained at 
large scales where as a disproportionately large fraction of the computational effort is 
expended on representing the small scales in fully-resolved simulations of such flows 
(e.g., see Pope, 2000-). Large eddy simulation (LES) is a technique that aims to 
explicitly capture the large, energy-containing scales while modeling the effects of 
the small scales that are more likely to be universal. This technique is both popular 
and is by far the most successful approach to modeling turbulent flows. We note, 
however, that in complex, wall-bounded and realistic configurations (such as, e.g., 
encountered in industrial situations), computational requirements for LES is still 
prohibitive that a hybrid (Reynolds Averaged Navier Stokes) RANS-LES approach 
is favored.- 

The nature of the dynamics of large scale circulation in the world oceans and 
planetary atmospheres is quasi two dimensional due to constraints of geometry (small 
vertical to horizontal aspect ratio), rotation and stable stratification. For example, 
consider the (inviscid and unforced) quasi-geostrophic equations that describe the 
dynamics of the large, geostrophically and hydrostatically balanced, scales: 



where q is potential vorticity approximated in the quasi-geostrophic limit by 



and u is the advection velocity approximated in the quasi-geostrophic limit by the 
geostrophic velocity given in terms of streamfunction %p hy \i ^ Ug = kx V-?/^. In other 
notation, V is the horizontal gradient operator, / is the Coriolis parameter given by 
f = fo + Py in the /3-plane approximation, y is the meridional coordinate, A^^ is 
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the Brunt- Vaisala frequency given in terms of the specified density gradient by A^^ = 
On the one hand, the particle-wise advection of potential- vorticity, and the dual 
conservation of (quadratic quantities) energy and potential-enstrophy are properties 
shared by quasi-geostrophic dynamics in common with two-dimensional flows. On the 
other hand, quasi-geostrophic dynamics shares in common with three-dimensional 
flow, the property of vortex stretching (in this limit, it is only the planetary vorticity 
/o that is stretched and is represented by the d/dz term in (j2]). 

It is the qualitative similarity of turbulence in these systems with two-dimensional 
turbulence, as elucidated by Charney, 1971,^ that is the primary reason for interest in 
two-dimensional turbulence. The dual conservation of (potential) enstrophy and en- 
ergy in (quasi) two-dimensional turbulence leads to profound differences as compared 
to fully three-dimensional turbulence: there exist two inertial regimes — a forward- 
cascade of (potential) enstrophy regime and an inverse- cascade of energy regime — in 
(quasi) two-dimensional turbulence in contrast to the single forward-cascade of en- 
ergy regime in fully three-dimensional turbulence. 

In the context of LES, which aims to model the effects of small-scales, it is clearly 
the forward- cascade inertial regimes that are of direct relevance. One of the most 
popular LES model is the Smagorinsky model,- and this class of eddy-viscosity mod- 
els assume that the main effect of the unresolved scales is to remove, from the resolved 
scales, either energy for 3D flows or (potential) enstrophy for (quasi-geostrophic) 2D 
flows — the appropriate quantity that is cascading forward. However, an examina- 
tion of the statistical distribution of the transfer of either energy in 3D turbulence 
or (potential) enstrophy in (quasi-geostrophic) 2D turbulence in the forward cas- 
cade regime^""- demonstrates that the net forward cascade results from the forward- 
scatter being only slightly greater than the backscatter. Clearly, models such as the 
Smagorinsky model, or more generally scalar eddy-viscosity, by modeling only the 
net forward cascade, fail to represent possible important dynamical consequences of 

3 



backscatter. 

The recent reinterpretation of the classic work of Leray — which considered a math- 
ematical regularization of the advective nonlinearity — in terms of the LES formalism, 
has given rise to the so-called regularization approach to modeling turbulence. ^'S-' 
An important model in this approach is the Lagrangian-Averaged Navier-Stokes-a 
(LANS-a) model introduced by Holm and co-workers.- 

The origins of the LANS-a turbulence model lie in 1) the notion of averaging 
over a fast turbulent spatial scale a, the reduced-Lagrangian that occurs in the 
Euler-Poincare formalism of ideal fluid dynamics^, and in 2) three-dimensional gen- 
eralizations^ of a nonhydrostatic shallow water equation system, known in literature 
as the Camassa-Holm equations.^- However, viewed from the point of view of the reg- 
ularization approach, this model can be thought of as a particular frame-indifferent 
(coordinate invariant) regularization of the Leray type that preserves other impor- 
tant properties of the Navier-Stokes equations such as having a Kelvin theorem. To 
add to the richness of this model, almost exactly the same equations arise in the 
description of second-grade fiuidsii^i^ and vortex-blob methods. ^'S-' 

There is now an extensive body of literature covering various aspects of the LANS- 
a model. In particular with respect to its turbulence modeling characteristics, analyt- 
ical computation of the model shear stress profiles has shown favorable comparisons 
against laboratory data of turbulent pipe and channel fiows^^ and a posteriori com- 
parisons of mixing in three-dimensional temporal mixing layer settings,^ in isotropic 
homogeneous turbulence settings,—*^ and in anisotropic settings^>^ compare well 
against Direct Numerical Simulations (DNS). In three dimensions, it has, however, 
recently been noted^^ that the use of LANS-a model as a subgrid model can be 
deficient in certain respects. In the two dimensional and quasi-two dimensional con- 
texts, a posteriori comparisons of LANS-a based computations have shown favorable 
comparisons against eddy-resolving computations.—""— Nevertheless, this model has 
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mostly been viewed as a complementary approach to modeling turbulence. 

The nonlinear gradient model^i"^ and the Rational LES model^~— are part of 
another class of LES models, built on a direct dynamical analysis of what should be 
a good approximation of the effect of the subgrid scales on the largest scales, through 
turbulent stresses. When the large scale velocity field is defined by low-pass filtering 
the velocity field, a natural asymptotic expansion leads to approximated turbulent 
stresses. This defines the nonlinear gradient model. An essential point is that the 
actual turbulent stresses of 2D and quasi-geostrophic turbulent flows, computed from 
direct numerical simulation, have been shown to be well approximated by the one 
defining the nonlinear gradient model.""^ 

The nonlinear gradient model flTT]) uses a natural approximation of the turbulent 
stresses. However this model has several drawbacks. Indeed, whereas it has been 
proven that the nonlinear gradient model turbulent stress (ITT!) preserves energy for 
two dimensional flows,— this is generally not the case in three dimensional flows, and 
instabilities or flnite time energy blow up can occur. The situation is not much better 
in 2D and quasi-geostrophic flows in that the incompressiblity constraint implies 
that the divergence of the deformation tensor {a in equation (fTTj) ) generally has a 
positive definite direction and a negative definite direction. Physically, this amount 
to an anistropic viscosity with positive value in some directions and negative values 
in other directions.-"-*^ These drawbacks mean that the nonlinear model is not 
a good physical model and will lead to instabilities, for two dimensional, quasi- 
geostrophic and three dimensional fiows. An alternative model based on entropic 
closures, keeping the main properties of the nonlinear model (good approximation 
of the turbulent stresses, conservation of energy), has been proposed and proven to 
give very good results for two-dimensional fiows.— In three dimensions, Domaradzki 
and Holm, 2001,— note that one component of the LANS-a (subfilter stress) model 
corresponds to the subfilter stress that would be obtained upon using an approximate 
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deconvolution procedure on the nonlinear gradient model. 

Analysis of the drawbacks of the nonlinear gradient model led Galdi and Layton 
to propose the Rational LES model. The Rational LES model coincides with the 
nonlinear model at leading order, but provides a stronger attenuation of the smallest 
scale. As confirmed by recent mathematical results,— the Rational LES model is well 
posed and should lead to stable numerical algorithms. It is thus a good candidate 
for LES. 

The Nonlinear- Gradient model has been well studied over more than three 
decades. These studies started with Leonard, 1974^ and Clark, Ferziger and 
Reynolds, 1979.^^ Rather than attempt an incomplete survey of the literature rel- 
evant to the a priori and a posteriori testing of this model here, we note that a 
fairly modern account of this can be found in Meneveau and Katz, 2000.— The 
more recent aspect of the Rational LES model is in making the highly favorable a 
priori comparisons of the Nonlinear- Gradient model more amenable to a posteriori 
simulations. For example, Iliescu et al., 2003^^ compare the behavior of the Rational 
LES model to the Nonlinear Gradient model (and the Smagorinski model) in the 
2D and 3D cavity flow settings, both at low and high Reynolds numbers. They 
find a) that laminar flows are correctly simulated by both models, and b) that at 
high Reynolds numbers, the Nonlinear Gradient model simulations, either with or 
without the Smagorinski model, lead to a finite time blow-up while the Rational 
LES model simulation displayed no such problem and succeeded in its LES role, i.e., 
compared to a fine-scale resolved simulation, the Rational LES model was able to 
capture and model the large-eddies well on a coarse mesh. Furthermore, they find 
that the Rational LES model performed better than the Smagorinski model alone 
in capturing the behavior of the large-eddies. Finally, we note that with both the 
Rational LES model and the LANS-a models, the burden of modeling borne by the 
additional dissipative term is smaller than in other approaches. 
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Following the development of these models, we note that the LANS-a and the 
rational LES model have interesting complementary properties: While the LANS-a 
preserves the Euler equation structure through the Kelvin theorem, the Rational 
LES model develops a good approximation of the turbulent stresses while ameliorat- 
ing problems associated with the nonlinear gradient model. It would thus be useful 
to examine the relation between these two models. In this article, we demonstrate 
the equivalence of the LANS-a to the Rational LES model in two-dimensions. By 
equivalence, we mean here that, the evolution equations for one of the models can 
be exactly transformed into the other. As will be evident, given the very differ- 
ent approaches taken in arriving at these models, it'll involve more than a simple 
transformation; it will also involve disentangling the turbulence term implied by the 
particular regularization of the nonlinear term. The importance of this result lies 
in the fact that mathematical results obtained for one of these models become also 
true for the other. We also demonstrate that these two models are different in three 
dimensions. 

In sections 1 and 2, after recalling the framework of turbulent stresses and LES, we 
briefly describe the nonlinear gradient, the Rational LES, and the LANS-a models. 
In section 3, we prove the equivalence of the Rational LES and of the LANS-a 
models in two dimensions. In section 4, we prove that they are not equivalent in 
three dimensions. After a brief numerical example, implications of the above results 
are discussed in the final section. 

II. LES OF TWO-DIMENSIONAL TURBULENCE AND THE 
NONLINEAR-GRADIENT MODEL 

In LES, the resolution of energy containing eddies that dominate flow dynamics is 
made computationally feasible by introducing a formal scale separation.- The scale 
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separation is achieved by applying a low-pass filter G with a characteristic scale a 
(2a^ is the second moment of G) to the original equations. To this end, let the 
fields u, g, etc... be split into large-scale (subscript Z)and small-scale (subscript s) 
components as 

U = + u^, 



where 



u;(x) = / G(x — x')u(x')(ix', 

J D 



U^Xj = U - Ui, 

the filter function G is normalized so that 



G(x')rfx' = 1, 

D 

and where the integrations are over the full domain D. In contrast to Reynolds 
decomposition, however, generally, uu ^ and u^; 7^ 0. 

For convenience, we write the two-dimensional vorticity equation as 
ID d 

-^ = ^ + u-Vu = V xF2d + V xD2d = F + D, (3) 
Ut ot 

where F2d is the two-dimensional momentum forcing, D2d is dissipation, and where, 
for brevity, we denote V x F^d by F, and V x D-^d by D. Applying the filter to ([3]) 
leads to an equation for the evolution of the large-scale component of vorticity which 
is the primary object of interest in LES: 

^ + V ■ {Miui) = Fi + Di-V-(T (4) 

OT 

where 

a = {uu)i - uiui (5) 

is the turbulent sub-filter vorticity-flux, and as in ([3]), we denote (V x F2d)i by Fi, 
and so also for dissipation. This turbulent subgrid vorticity-flux may in turn be 
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written in terms of the Leonard stress, cross-stress, and Reynolds stress^ as 

a = {uiui)i - uiui + {uiUs)i + {usUi)i + {usUs)i . (6) 



Leonard stress Cross-stress Reynolds stress 

However, while a itself is Galilean-invariant, the above Leonard- and Cross-stresses 
are not Galilean-invariant. Thus when these component stresses are considered indi- 
vidually, the following decomposition, originally due to Germano, 1986^^ is preferable 

a = {uiuji)i - uiiuu + {uiujs)i + insUJi)i - UaUsi - UsiUu + (u^w^); - ^siUJsi ■ (7) 

^ V ' ^ v ' ^ V ' 

Leonard Stress Cross-stress Reynolds stress 

The filtered equations, which are the object of simulation on a grid with a res- 
olution commensurate with the filter scale in LES, are then closed by modeling 
subgrid-scale (SGS) stresses to account for the effect of the unresolved small-scale 
eddies. In this case (jlj) will be closed on modeling the turbulent subgrid vorticity-fiux 
a. 

As is tradition, a Gaussian filter is chosen. In eddy-permitting simulations, some of 
the range of scales of turbulence are explicitly resolved. Therefore, information about 
the structure of turbulence at these scales is readily available. In LES formalism, 
there is a class of models that attempt to model the smaller unresolved scales of 
turbulence based on the assumption that the structure of the turbulent velocity field 
at scales below the filter scale is the same as the structure of the turbulent velocity 
field at scales just above the filter scale.— 

Further expansion of the velocity field in a Taylor series and performing filtering 
analytically results in 

a quadratic nonlinear combination of resolved gradients for the subgrid model.— 
The interested reader is referred to Meneveau and Katz, 200022. for a comprehensive 
review of the nonlinear-gradient model. 
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Equivalently, expansion of u/ and ui in the Galilean-invariant form of the Leonard- 
stress component of the sub-filter eddy-flux of vorticity ([7]) in a Taylor series: 



,UiUJi I - UiiUii 



= I G(x - x') (^u,(x) + {x' - ^)i|^(x)) [u^i{^) + {x' - ^).|^ - 

j dK' G(x-x') (^uKx) + {x' - x),|^(x)^* j rfx' G(x-x') l^uii^) + {x' - a:),|^(x)) 
produces at the first order 

a = 2a'p^ ^ + 0{a') = 2a'Vui ■ Vui + 0{a'), (9) 

OXj OXi 

where is the second moment of the filter used. The leading order is again a 
quadratic nonlinear combination of resolved gradients. The approximate model that 
retains only the second order term is called the nonlinear gradient model. In this 
two-dimensional setup, it reads 

^ + ui-Vuji = -2c? [Vuf .V(Va;0] + i^z + A (10) 

(please see the appendix for the definition of operator Vuf.V). 

For simplicity, we have presented the two-dimensional derivation of the nonlin- 
ear gradient model, however similar considerations lead to the three dimensional 
nonlinear gradient model: 

^ + . Vui = -2^2 V • [Vui Vui] - VP + {F^d\ + [D^i), . (11) 

In the two dimensional context, this model has been derived by Eyink, 2001— without 
the self-similarity assumption, but rather by assuming scale-locality of contributions 
to a at scales smaller than the filter scale, and its use has been investigated by various 
authors.-i^ Nadiga, 2008 and 2009^'^ have demonstrated excellent a priori testing 
of the nonlinear gradient model in quasi-geostrophic turbulence, the same also holds 
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in the three-dimensional turbulence context. ^'S-' ^ The nonlinear gradient model, 
however, holds much better in two-dimensional and quasi two-dimensional settings 
than in fully three-dimensional settings. 



III. RATIONAL LES MODEL AND THE LANS-a MODEL 
A. Rational LES model 

By analyzing the nonlinear-gradient model in terms of Fourier components, Galdi 
and Layton, noted that the nonlinear-gradient model increases the high wavenumber 
components (scales smaller than the filter scale) and therefore does not ensure that 
Ui is smoother than u. Consequently, to remedy this problem, they proposed an 
approximation which attenuates the small scale eddies, but is of the same order 
accuracy for large eddies (the two approximations coincide at order a^, see ( ITTj) ). 

To this end, rather than using a Taylor expansion of the filter [e~^^^ 1 — bx^), 
they considered the rational approximation 



Using the above sub-diagonal Fade approximation, the modified nonlinear-gradient 
model leads to the 'Rational LES' model. We refer to Galdi and Layton, 2000^^ for 
the derivation of the evolution equation for u; (which is an approximation of the 
large scale component of the full velocity field u.) It is 



^ + • Vu; = -2a' (/ - a' A) ' V • [Vu^Vu^] - VP + (F3,), + {D^d), , (13) 




1 




1 + 6x2 



with V • = 0, and where (/ — a'' A) is the inverse of the operator (/ — a 
(easily expressed in a Fourier basis). 



:2A) 
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B. The LANS-a model 



In the context of the three-dimensional incompressible Navier-Stokes equations 

— + u • Vu = -V0 + F^d + «^Au ; V ■ u = 0, (14) 
at 

on a suitable domain with appropriate boundary conditions, Leray regularization of 
is expressed by ^'S-' ^ 

dxi 

— + ■ Vu = -V0 + Fsrf + uAu, (15) 

CM 

where is the large scale component of velocity filtered at a characteristic length a, 
(f) = p/po is the normalized pressure, is the external forcing and u the kinematic 
viscosity. The filtered velocity can be obtained by application of a convolution 
filter to u. A particularly important example is the Helmholtz filter, to which we 
turn momentarily. The Leray approach is basic to many recent studies in regularized 
turbulence. This regularization model does not preserve some of the properties of 
the original equations (fT^ . such as a Kelvin circulation theorem. This is where the 
LANS-a formulation provides an important extension. 

A transparent way to present the LANS-a model is obtained when the incom- 
pressible Navier-Stokes (momentum) equations are written in the equivalent form 

^ - u X (V X u) = -V0 + + i^Au. (16) 

The LANS-a model is then given by ^'S'' -"^ 

^ - X (V X u) = -V0 + + lyAu. (17) 

Thus, just as the Leray regularization corresponds to the filtering of the advecting 
velocity, the LANS-a regularization amounts to filtering the velocity in the nonlinear 
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term when written as the direct product of a velocity and a vorticity lj = V x u. 
The LANS-a model may be written in the more common advective nonlinearity form 

^ + Ui ■ Vu - {Vuif Aui = -Vp + F^a + ^^Au. (18) 

The filtered velocity is obtained by an inversion of the Helmholtz operator: u/= (1 — a^A)^^u 
with appropriate boundary conditions. (It has to be noted that in a non-periodic 
domain, the boundary conditions that are necessary to invert the Helmholtz oper- 
ator are specific to this modeling procedure.) The third term on the left in ( !T8|) is 
introduced in the LANS-a modeling approach to restore a Kelvin theorem to the 
modeled equations. 

It is also instructive to consider the evolution of vorticity. For the Navier-Stokes 
equation, vorticity evolution is 

— + u • VcJ = a; ■ Vu + V X + z/Au;. (19) 

The vorticity evolution corresponding to the LANS-a model is 
d(jj 

— + ui-VoJ = u ■ Vui + V X Fsrf + uAuj, (20) 
ot 

where in addition to a filtered advecting velocity, a mollification of the vortex- 
stretching term is evident. 

In two dimensions, (120!) reduces to 

^ + ui-Vuj = V X F2d + V X D2d = F + D, (21) 



where forcing and dissipation have been written in correspondance with the notation 
used in the two dimensional vorticity equation ([3]) and its LES counterpart (jlj). 
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IV. IDENTITY OF THE RATIONAL LES AND LANS-a MODELS 
IN TWO DIMENSIONS 



In this section, we consider the Rational LES model and the LANS-a models in 
two- dimensions. Taking the curl of the two-dimensional velocity equation for the 
Rational LES model ( IT3l) . we obtain the vorticity equation 

^ + ui-ViOi = -2a^ (/ - a^Ay' [Vuf .V{VuJi)] + Fi + (22) 

where ui is the vertical component of uJi. In order to compare the Rational LES 
model ([22]) with the LANS-a model we apply operator (J - a^A) to and 
write the evolution equation for u as 

^ + Vii -Vuj = 6M + F + D. (23) 
at 

Comparing fl23|) with fl2Tl) . we note that 5M is the difference between the two models 
and is given by 

5M = -2a^ [Vuf.V{Vuji)] + • V [(/ - a^A) ui] - (/ - a^A) [m • Vooi] . 

By direct computation this expression simplifies to 

6M = o?{~2 [Vuf .V(Vu;0] - u • V \Aui\ + A [u^ • Vt^^]} . 

Then using the vector calculus identity fl2^ in the appendix, we conclude that bM = 
0. The dynamics of u is thus the same as given by the LANS — a model 

did ^ ^ 

— + ui-Vu = F + D. 

ot 

We thus conclude that the Rational LES model and the LANL-a models are equiv- 
alent in two dimensions. 
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V. THE RATIONAL LES AND LANS-a MODELS ARE DIFFERENT 
IN THREE DIMENSIONS 



In three dimensions, the Rational LES model for an incompressible flow (V.u = 0) 

is 

^ - X (V X uO = - 2a^ (/ - a^A)"' [Vuf .V(VuO] + + u/\mu 

^ (24) 
where Pi is the sum of the physical and kinetic pressure. The LANS-a model is 
(9ii VP 

— - X (V X u) = - + F^, + uAu, (25) 

at p 

where u = (/ — a^A) u;. 

Applying the operator (/ — a^A) to (!24l) . we obtain the equation verifled by u = 
(/ — a^A) u; in the case of the Rational LES model: 
du VP^ 

— - X (V X u) = ^ + N + + uAu, (26) 

ot p 

with = {I - a^A) P2, and with 

N = {-2Vuf .V(Vui) + u; x V x [Au^] + A [u; x (V x u^)]} . 

The two equations (|25ll and ( !26l) are equivalent if and only if N is a gradient, that 
is if and only if V x N = 0. For two-dimensional vector-flelds, we have proven in 
the section above that 5Me^ = V x N = 0. In contrast this is wrong in general for 
three-dimensional vector flelds, because of vortex-stretching type terms present for 
three dimensional vector flelds and non-present for two dimensional vector flelds. In 
order to prove this we give below an example of fleld u; for which V x N 7^ 0. 

Consider for example u; = y'^ex—xzey+xyez. Then ur^ is actually non- divergent: 
V ■ u; = 0. By direct computation, we have V x [Au;] = 0, u; x V x [Au^] = 0, 
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V X [Vuf.V(Vui)] = 2e^, V x {A [u; x (V x u/)]} = -8e^ and then V x N = 

4a2e, ^ 0. 

We thus conclude that the Rational LES and the LANS-a models are not equiv- 
alent in three dimensions. 

VI. A NUMERICAL EXAMPLE 

The primary emphasis of this article is the above analytical demonstration of the 
equivalence of the Rational LES model and the LANS-« model in two dimensions 
rather than an evaluation of the performance of the model(s) considered. Neverthe- 
less, at the insistence of one of the referees, we briefly present an example computa- 
tion in two dimensions in this section. 

We consider a (stochastically) forced-dissipative flow in a doubly period domain 27r 
on the side. As is conventional in numerous earlier investigations of two dimensional 
turbulence, dissipation D consists of linear damping: — rcj, where r(= 10^^) is a 
frictional constant, and an eigth order hyperviscous term that acts as a sink of the 
net-forward cascading enstrophy. Forcing F is scaled as F = \/2tF^ where F is an 
isotropic stochastic forcing in a small band of wavenumbers 15 < /c/ < 16 drawn 
from an independent unit variance Gaussian distribution and which is temporally 
uncorrelated: (^Fy_(t)Fk/(t')^ = 5kk'(i — t')- A fuUy-deahased pseudo spectral spatial 
discretization is used in conjunction with an adaptive fourth-fifth order Runge-Kutta 
Cash-Karp time stepping scheme. The time step used ensures that the relative error 
of the time increment is less that 10~^, and with the time step ending up being much 
smaller than that required by stability requirements. 

For the reference computation, a 512x512 physical grid is chosen giving a gridsize 
of 7r/256. Figure 1 shows the vorticity field after the fiow has equihbrated (at t=2600 
eddy turnover times.) Given the stochastic forcing and the turbulent nature of the 
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100 200 300 400 500 

x-coordinate 

FIG. 1. Snaphot of vorticity field of reference run at statistical stationarity. The stochastic 
forcing is confined to a band of wavenumbers between 15 and 16 (domain is 2tt x 2tt) and 
dissipation consists of a combination of Rayleigh friction and hyperviscosity. 

flow, a statistical consideration of the flow is in order: The evolution of the domain 
integrated kinetic energy and enstrophy as a function of time is shown in Fig. 2. 
Figure 3 shows the one-dimensional spectral density of kinetic energy and the spectral 
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FIG. 2. The evolution of kinetic energy and enstrophy (inset) in the reference run (sohd 
hne; black), bare truncation run (dashed line; red) and the LES run (dot-dashed line; 
green). The horizontal axis in the inset spans the same range of times as in the main 
plot. The former uses a 512x512 physical space grid where as the latter two runs use a 
128x128 grid. The changes introduced by the subgrid model are so as to improve the bare 
truncation run in the direction of the reference run. 

flux of kinetic energy u*-u-Vu) where Re(.) denotes the real part,? 

denotes the Fourier transform and superscript * denotes complex conjugate.) For 
the spectral flux of kinetic energy diagnostic, we verify that the integral over all 
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wavenumbers goes to zero to within roundoff. Note that the x-axis in Fig. 3 is 
truncated at wavenumber 60 to better focus on the range of scales of interest. In all 
these figures, the solid line represents the reference run. 

Next, we choose a filter width of 7r/32, and following arguments similar to those 
in section 13.2 of Pope, 2000^, we choose an LES gridsize of 7r/64 (that corresponds 
to a 128x128 physical grid.) On the (coarser) LES grid, we perform two simulations: 
One that we call a bare truncation — fl22|) . but without the first term on the right 
hand side-and a second one with the LES model discussed above — fl22l) — with the 
rest of the setup being identical. The evolution of domain-integrated kinetic energy 
and enstrophy and the spectral density and spectral flux density of kinetic energy for 
these two simulations are shown again in Figs. 2 and 3. The bare truncation run in 
these figures is shown by dashed (red) lines, whereas the Rational LES or LANS-a 
model runs are shown by dot-dashed (green) lines. 

In each of these diagnostics, the tendency of the model to improve on the bare 
truncation is evident. In the spectral density plot, on comparing the bare truncation 
with the model simulation, the tendency of the model to deemphasize the small 
scales while increasing the energy in the large scales is seen. The dynamics of how 
this is achieved is seen to be that of an augmentation of the inverse cascade by the 
model term as indicated by the blue line in the spectral flux inset. The net result 
is that the full nonlinear flux of energy shows a smaller forward cascade^^ and an 
increased inverse cascade, as compared to the bare truncation simulation. And these 
changes in the spectral flux of energy are in the direction of bridging the (coarser) 
bare truncation run to the reference simulation. We note that a) we did not tune any 
of the parameters to match the reference run; we anticipate that with tuning, the 
LES results could better match the reference run, and that b) the computations on 
the LES grid are about 60 to 100 times less computationally intensive as compared 
to the reference run, with the overhead for the model (over bare truncation) being 
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One dimn. wavenumber 

FIG. 3. The one-dimensional power spectral density (logarithmic scale) and spectral flux 
density (inset; linear scale) plotted as a function of the one-dimensional wavenumber (log- 
arithmic scale). The horizontal axis in the inset spans the same range of wavenumbers as 
in the main plot. Reference run: solid line, black; bare truncation run: dashed line, red; 
and the LES run: dot-dashed line, green. The thick dot-dashed line in blue in the inset 
corresponds to the spectral flux of energy due to the subgrid model. In the range of scales 
where there is an inverse cascade of energy, the LES run is more energetic than the bare 
truncation run, and the LES run closely follows the reference run. However at the small 
scales, the energy spectrum of the LES run falls off faster than the bare truncation run (at 
these range of scales, the reference run is still inertial.) The increased level of energy at 

the large scales in the LES run is seen as due to a secondary inverse cascade that is put 
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in place by the subgrid model (backscatter). In effect, as compared to the bare truncation 
run, in the LES run, the forward cascade of energy is reduced and the inverse cascade of 
energy is augmented. 



negligible. 



VII. CONCLUSION 

In its popular form, the LES approach to modeling turbulence comprises of ap- 
plying a filter to the original set of equations; the nonlinear terms then give rise 
to unclosed residual terms that are then modeled. However, the regularization ap- 
proach to modeling turbulence consists of, besides other possible considerations, a 
modification of the nonlinear term based on filtering of one of the fields. The lat- 
ter approach, however, implies a model of the unclosed residual terms when viewed 
from the point of view of the former. We consider the Rational LES model^ that 
falls under the former approach, and the LANS-a model^ that falls under the latter 
approach. In this article, we demonstrate that the two models are equivalent in two 
dimensions, but not in three dimensions. Their equivalence in two dimensions allows 
arguments about the mathematical structure and physical phenomenology of either 
of the models to be equally valid for the other. 
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VIII. APPENDIX 

We derive in this appendix calculus identities for two-dimensional or three- 
dimensional vector fields. We define for any vector fields A and scalar 5: VA^.V(V-B) = 
diAjdijB (sum over repeated indices, here and in the following) and AA-^.V-B = 

21 



AAjdjB. Similarly, for any two vectors fields A and B, we define VA"^.V(VB) = 
diAjdijB and AA^.VB = AAjdjB 

1. We first note that the useful, and easily derived, vector calculus identity 

A (A.V5) = AA^.VB + 2VA^.V(Vfi) + A.V (AB) . (27) 

2. We then note that for a 2D solenoidal vector field u (V.u = 0), if a; = 
(V X u) .e^, then 

Au^.Vw = 0. (28) 

Indeed using the flow incompressibility V.u = 0, we have dxoo = Auy and 
dyU = —Aux, then Au-^.Vo; = AuxdxOJ + AuxdxOJ = 0. 

3. Then for a 2D incompressible vector field u with w = (V x u) .e^, using f l27p 
and ( l28l) we obtain 

A (u.Vw) = 2Vu^.V(Va;) + u.V (Aw) . (29) 
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1 Snaphot of vorticity field of reference run at statistical stationarity. 
The stochastic forcing is confined to a band of wavenumbers between 
15 and 16 (domain is 2tt x 2tt) and dissipation consists of a combination 
of Rayleigh friction and hyperviscosity 



2 The evolution of kinetic energy and enstrophy (inset) in the reference 
run (solid line; black), bare truncation run (dashed line; red) and 
the LES run (dot-dashed line; green). The horizontal axis in the 
inset spans the same range of times as in the main plot. The former 
uses a 512x512 physical space grid where as the latter two runs use a 
128x128 grid. The changes introduced by the subgrid model are so as 
to improve the bare truncation run in the direction of the reference run 
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The one-dimensional power spectral density (logarithmic scale) and 
spectral flux density (inset; linear scale) plotted as a function of the 
one-dimensional wavenumber (logarithmic scale). The horizontal axis 
in the inset spans the same range of wavenumbers as in the main 
plot. Reference run: solid line, black; bare truncation run: dashed 
line, red; and the LES run: dot-dashed line, green. The thick dot- 
dashed line in blue in the inset corresponds to the spectral flux of 
energy due to the subgrid model. In the range of scales where there 
is an inverse cascade of energy, the LES run is more energetic than 
the bare truncation run, and the LES run closely follows the reference 
run. However at the small scales, the energy spectrum of the LES run 
falls off faster than the bare truncation run (at these range of scales, 
the reference run is still inertial.) The increased level of energy at 
the large scales in the LES run is seen as due to a secondary inverse 
cascade that is put in place by the subgrid model (backscatter) . In 
effect, as compared to the bare truncation run, in the LES run, the 
forward cascade of energy is reduced and the inverse cascade of energy 
is augmented 
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